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II superconductor picture, and with an effective dual gluon mass of 0.905 ± 0.163 GeV. 
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Figure 1: (left) Simple Wilson loops that make the gqq Wilson loop, (right) Our U shape gluon-quark- 
antiquark geometry. 

1. Introduction 

Here, we compute in Lattice QCD the colour fields produced by static colour sources, with 
the aim to study the dual superconductor analogy of QCD. Detailed studies of the field profile of 
the tubes, including the string merging/separation, the Casimir scaling, and the penetration in the 
vacuum, are presented. 

In section 2, we review the calculation of chromo-fields in the lattice, using the Wilson Loop. 
Then, in section 3 the numerical results are given for different cases. Namely, for flux tube profile of 
the static gqq in a U shape geometry. For the Casimir coefficientes of the different representations 
of SU (3) and for the dual superconductor parameters of QCD. We conclude in section 4. 

2. The Wilson Loops and Colour Fields 

The Wilson loop for the static hybrid gqq was deducted in [1, 2, 3, 4], therefore we only present 
the fundamental expressions. The Wilson loop for this system is given by 

W gq q = W 1 W 2 -^W 3 (2.1) 

where W\, W2 and W3 are the simple Wilson loops shown in Fig. 1. 

We also study the interaction between two charges of different representations of SU (3) in a 
colour singlet. So for each different SU (3) representation we will have a different Wilson Loop. 



The different Wilson loops are [5], 

W 3 = trU (2.2) 

W 8 = |Ws| 2 - 1 (2.3) 

W 6 = ^((trC/) 2 + trC/ 2 ) (2.4) 

W l5a = trUW 6 -trU (2.5) 

Wio = ^ ((trC/) 3 + 3trC/(trC/) 2 + 2trC/ 3 ) (2.6) 

W 2 4 = trUW 10 -W 6 (2.7) 

W21 = \W 6 \ 2 -\W 3 \ 2 (2.8) 

W 15s = — ((trC/) 4 + 6(trC/) 2 trC/ 2 + 3(trC/ 2 ) 2 + 8trC/trC/ 3 + 6trC/ 4 ) (2.9) 
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Figure 2: Results for the U geometry at y = 4 and z = 0. 

In order to improve the signal to noise ratio of the Wilson loop, we use APE smearing, [4], 
with w = 0.2 and iterate this procedure 25 times in the spatial direction. To achieve better accuracy 
in the flux tube, we apply the hypercubic blocking (HYP) in time direction, [6], with a\ = 0.75, 
a 2 0.6 and a 3 0.3. 

The chromoelectric and chromomagnetic fields are given by, 

(zf) = (Poi)-^ (2-10) 

where the jk indices of the plaquette, P, complement the index / of the chromomagnetic field. 
Notice that we only apply APE and HYP to the Wilson loop and not to the plaquette P. 
The energy (^f 7 ) and lagrangian densities are given by 

±((E 2 ) + (B 2 )) (2.12) 

JSf = \ ((E 2 ) - (B 2 )) . (2.13) 

We will also try to fit the results to the solutions of the Ginzburg Landau equations of super- 
conductivity, which in Lorentz-Heaviside units are given by 

— (V-i^-AW + cW + bm 2 ^ = (2.14) 
2m nc 

VxVxA= ^^(^*V^-W^*-^|^| 2 A). (2.15) 
lime nc 
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Figure 3: On the left, results for r = y at x = and z = 0. On the right, for r = Vx 2 +z 2 at the mediator 
plane. 

This equations define two length scales, the coherence length t; and the penetration length A, 




(2.16) 



3. Results 

Here we present the results of our simulations with 286 24 3 x 48, j8 = 6.2 quenched con- 
figurations generated with the version 6 of the MILC code [7], via a combination of Cabbibo- 
Mariani and overrelaxed updates. The results are presented in lattice spacing units of a, with 
a = 0.07261 (85) fm or a' 1 = 2718 ± 32MeV. 

In this work a geometries for the hybrid system, gqq, is investigated: a U shaped, defined 
in Fig. 1. We also investigate the case of two static colour charges only, but for different SU(3) 
representations, staring by the wuark-antiquark and gluon-gluon system, but we also study higher 
representations up to the 15c . 

The use of the APE (in space) and HYP (in time) smearing allows us to have better results for 
the flux tube, while suppressing the fields near the sources. 

3.1 U profiles 

In Fig. 2 we present the profiles for the U geometry for I = 8 and d between and 16 at 
y = 4. We can see the stretching and partial splitting of the flux tube in the equatorial plane (y = 4) 
between the quark and antiquark. For d < 4 the separation between the two flux tube is not visible. 
This is due to the overlap of the tails of the flux tube profile, i. e. to the width of the flux tubes. 
For larger separations the two flux tubes maxima are separated, while the tails of the flux tubes 
contribute to a non zero field at x = 0. 

3.2 Casimir scaling 

We measure the quotient between the energy density of the two gluon glueball system and of 
the meson system, in the mediator plane between the two particles (x = 0). The results are shown 
in Fig. 3. 
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Figure 4: Flux tube profiles for different representations of SU(3). On the left, results for x = 0. On the 
right, for the mediator plane, with r = \fx 2 +z 2 . 







r~i~ s 

o 15a 
^ 10 

o 27 
"?~ 24 
15s 


jrr~r=f # ^ 


i- ; 

, 1_ 


. o #= o — -o j 

h-i — b^4^ — E=J 







6 s::: 



i k 



8 o 





8 


~$~ 


6 


o 


15a 




10 


o 


27 




24 


1 ~$~ 


15s 



y, x=0, z=0 



0.5 1 1.5 2 2.5 3 3.5 4 

r,y=0 



Figure 5: Ratios between the static potential of superior SU(3) representations to the fundamental one. The 
results are compared with the Casimir Scaling predictions ( black interrupted horizontal lines ). On the left 
the results are for x = 0. On the right they are for the mediator plane. 



We make a constant fit to the data in Fig. 3, the result for r = y at x = z = is 2.25 ± 0.02 
and for r = (x,z) at y = 4 is 2.24 ± 0.06. As can be seen, these results are consistent with Casimir 
scaling, with a factor of 9/4 between the energy density in the glueball and in the meson. This 
corresponds to the formation of an adjoint string. 

In Fig. 4 we present the flux tube profiles of the energy density for different representations 
of SU (3). In this case, we apply three levels of HYP in the temporal direction, in order to further 
reduce the signal to noise ratio. 

In Fig. 5 we see the ratios between the results for the flux tube profiles in the different higher 
representations of SU (3) and the fundamental one. The results presented are for short distancies, 
and seem to agree with Casimir Scaling predictions. However, for larger distances the statistical 
errors, become greater impeding us from concluding about the validity of the Casimir scaling 
hypothesis for the flux tube profiles. 



3.3 Dual Gluon Mass 

In 1970's, Nambu [8], 't Hooft [9] and Mandelstam [10] proposed a very interesting idea: 
quark confinement would be physically interpreted using the dual version of the superconductivity, 
where the QCD vacuum state would behave like a magnetic superconductor. The chromoelectric 
field originated by a qq pair would be squeezed by a Meissner effect into a dual Abrikosov flux tube, 
or even to QCD strings in the strong coupling limit, giving rise to the confining linear potential. 
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Figure 6: Fits of the profiles of the flux tubes of a meson (left) and of a glueball (right). 



Color confinement could be understood as the dual Meissner effect. 

Notice that in common superconductivity the magnetic field decays with B ~ e~ r ^ L and this 
can be interpreted in terms of an effective mass for the photon m r = I/Xl. There is also evidence 
for the dual superconductor picture from numerical simulations or theoretical modles of QCD 
[12,11,13]. 

Here we test two functions, ae~ 2 ^ r and clKq (jiir) where Ko is the modified Bessel function 
of order zero, and where \i = j-, Xl is the penetration length. So, in this case we have ji as a 
possible definition mass for the dual gluon. Fitting the chromoelectric field and the lagrangian 
density section in the mid distance of the flux tube of the meson and the two gluon glueball, we 
obtain the results presented in Table 1 for the effective dual gluon mass, of the order of ^ 1 GeV . 
Some values found in literature, for the effective dual gluon mass, [14, 15, 16, 17], and the effective 
gluon mass, [18], point to a similar value. Our fits are depicted in Fig. 6 and detailed in Table 1. 

4. Conclusions 

Notice that in type Type-II superconductors the flux tubes repel each other while in Type-I 
superconductors they attract each other and tend to fuse in excited vortices. In QCD this raises the 
question, how the two pictures, of one adjoint string and of two fundamental strings, with different 
total string tensions, match? 

By separating the quark from the antiquark in a static hybrid, we observe the fundamental 
string fusion and separation. When the quark and the anti-quark are superposed, we observe to the 
formation of an adjoint string between the two gluon and agrees with Casimir Scaling measured 
by Bali [5], This can be interpreted with a type-II superconductor analogy for the confinement in 
QCD with repulsion of the fundamental strings and with the string tension of the first topological 
excitation of the string (the adjoint string) larger than the double of the fundamental string tension. 

We further study the flux tube created by higher representations of SU(3), but above the de- 
cuplet 10 representation, it is not clear yet whether Casimir scaling continues to occur in the field 
profiles of the flux tube. The detailed study of the high representation field profiles still will con- 
tinue to require further and more detailed investigation. 

We also compute a value for the dual gluon mass estimated form the penetration length of 
the flux tube, to be X~ l = 0.91 ±0.16. It is remarkable that this result is computed with gauge 
independent lattice QCD calculations. 
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aK$(nr) 


M (GeV) 


X 2 /dof 


11 (GeV) 


X 2 /dof 


\ (r) 
(i)(«) K J 


1.170 ±0.228 


1.069 


0.805 ±0.287 


1.827 


^(\){a) (r) 


1.170±0.119 


0.512 


0.865 ±0.188 


1.203 


E (2)(a) W 


1.231 ±0.286 


1.547 


0.881 ±0.334 


2.084 




1.210±0.056 


0.887 


0.897 ±0.085 


1.185 




1.208 ±0.068 


0.560 


0.909 ±0.099 


0.909 


^(2)W M 


1.210±0.063 


1.162 


0.889 ±0.097 


1.262 


^(2)W (0 


1.191 ±0.031 


1.066 


0.899 ±0.048 


1.106 



Table 1: Results for the dual gluon mass, where (1) is for the two gluon glueball and (2) for the quark- 
antiquark cases, and (a) at y = 4 and z = with r = x and (b) at y = 4 with r = (jt, z) . 
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